clear all
count =1;
n=1;

for ii=1:20
 for jj=1:20
  A(ii,jj) = exp((ii-1)/4)*5e-5;
  AA=A(ii,jj)
  width(ii,jj) = exp((jj-1)/4)*1d-15;
  WW = width(ii,jj)
save('infile','n','AA','WW','-ascii')
system('a.out');
load('fort.111');
ent(ii,jj) = fort(end,6);
pop(ii,jj,1:3) = fort(end,2:4);
 end
end

